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ABSTRACT 

Ionized reflection has often been considered as the explanation for the unusual Fe Ka variabil- 
ity observed in MCG-6-30-15. In this paper, we test this model using a 325 ks observation of 
MCG-6-30-15 by XMM-Newton and BeppoSAX. The data are fit between 2.5 and 80 keV with 
the constant density models of Ross & Fabian. The best fit ionized reflection model requires 
the Fe Ka line to be split into two reprocessing events: one from the inner disc to build up the 
red wing, and the other from the outer accretion disc to fit the blue horn. The implied geom- 
etry is a disc which becomes strongly warped or flared at large radii. A good fit was obtained 
with a solar abundance of iron and a reflection fraction (R) of unity for the inner reflector. 
The combination of the two reflection spectra can appear to have R > 2 as required by the 
BeppoSAX data. The inner reflector has an ionization parameter log £ = 3.8, but the outer one 
is neutral with an inner radius ~ 70 gravitational radii (r fl ), corresponding to a light crossing 
time of about an hour for a 10 7 M Q black hole. Applying this model to time-resolved spectra 
shows that the inner reflector becomes more ionized as the source brightens. This reduces the 
strength of the red wing at high flux states. The X-ray source is constrained to arise from a 
narrow annulus at ~ 5 r g , with only 6 per cent of the 2-10 keV flux due to the outer repro- 
cessor. This amount of localized energy generation is extremely difficult to produce without 
resorting to other energy sources such as the black hole spin. In fact, all the Fe Ka models 
fit to XMM-Newton spectra of MCG-6-30-15 require a large increase in energy production at 
the inner edge of the accretion disc. 
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1 INTRODUCTION 

X-ray observations of active galactic nuclei (AGN) allow for an 
unparalleled probe of the physics at the inner edge of the accre- 
tion di^c J _Srjectr^_featoessuch as a strong Fe Ka line at 6.4 keV 
(see Reynol ds & Nowakl 120031) for a recent r eview) and a flatten- 
ing of the power-law slope above ~10 keV iPoundsetal. 1990; 
Nandr a& Pounds! 1994 are indications that some fraction of the X- 
radiation is bein g reprocessed (or reflected) in op tically-thick ma- 
terial (e.g jGeorge & Fabianll99lllMatt Perola & Piroll99lh . Us- 
ing the ASCA observatorv lTanakaet^lT^^95i) discovered a broad, 
asymmetric Fe Ka line in the bright Seyfert 1 galaxy MCG-6- 
30-15 (z — 0.008) that was consistent with the profile expected 
from an emission line produced n ear 6 gravitational radii (6 r Q ; 
r g = GM /c 2 ) from the black hole iFabian et alJl989HLaorll99lf) . 
Therefore, the reflection signatures observed in X-ray spectra of 
AGN are likely caused by reprocessing within the accretion disc - a 
plausible scenario if th e X-ray source resides in a hot co rona above 
the di sc <Galeev. Rosner & VaianJl979UHaardt & MaraschJl99lL 
Il993l) . 
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In the disc-corona geometry it is expected that the Fe Ka line 
will respond rapidly (~ 50r Q s for a 10 7 Mp black hole) to any 
change in the continuum <Stellall990LlMatt & PerobJl992h . Thus, 
the expectation is to observe a correlation between the line flux and 
the continuum flux, and a constant Fe Ka equivalent width (EW). 
However, many long observations of MC G-6-30-15 have shown 
that the Fe line exhibits little variability <Lee et all2000tlRevnoldsl 
l200d IShih. Iwasawa & Fabiadl2002l: iFabian et al]l2002bl) despite 
large (factors of 2-3) changes in the continuum flux, and an EW 
that decreases as the source brightens. Moreover, when variations 
were detected the y were found to be uncorrelated with the con- 
tinuum variability fawasawa et alll99dll999llVaughan & Edelsonl 
2001 ). This perceived lack of Fe Ka reverberation is a major chal- 
lenge to the disc-corona model of X-ray production in AGN. 

A possible solution to this problem is if the Fe l ine arises from 
the surface of a strongly photoionized accretion disc jlwasaw a et alj 
Il99d iLee et aljEoOO: iRevnoldslEbobl iBallantvne & Rossll2002h . 
When iron is highly ionized then increasing the strength of the il- 
lumination will result in a wea ker Fe Ka line re lative to the con- 
tinuum (see Fig. 2 in lMatt. Fabian & Rosslll996h . The flux of the 
line can thus remain approximately constant despite variations in 
the driving continuum. This situation can occur when helium-like 
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Fe is the dominant ionization state <Matt. Fabian & Rossi [l993t 
iBallantvne. Fabian & R oss 2002), and thus requires the disc sur- 
face to have a high ionization parameter, log f > 3, where £ = 
47rJ 7 x/nH- The major difficulty with this solution is that the He- 
like Fe Kq line occurs at 6.7 keV, but the line core of MCG-6-30- 
15 is well established t o be at 6.4 keV, and does not arise from very 
distant material iLee et all20 02). Therefore, at first glance it seems 
unlikely that the Fe Ka line arises from an ionized disc. 

One way to check this interpretation is to compare the 
entire spectrum to models of X-ray reflection from accre- 
tion discs. Calculations of X-ray reprocessing have shown 
that the shape of th e spectrum is sens i tive to the ionization 
state of the material iRoss & F abian||l993|:|Ross 1 Fabian & Yound 
1999), its density distribution jN ayakshin. Kazanas & Kallman 
20001 IBallantvne. Ross & Fabianl l200ll: iDumont et alj 1200? 
R6za hska^t^iT2002l). and the properties of the illuminating radia- 
tion iBallantvne & Ross 2002). These models self-consistently cal- 
culate the strength of the Fe Ka line which changes as a function of 
ionization state and abundance l Balla ntyne et all2 002). so the line 
and the entire continuum can be fit simultaneously with only a few 
parameters. 

T his exercis e was prev iously performed for MCG-6-30-15 by 
IBallantvne & Fabianl J200ll) using ASCA data. They found that the 
data between 3 and 10 keV (the 3 keV cut-off was used to avoid the 
warm absorber) was well fit by an ionized reflector, but required 
the addition of a separate line component further out along the disc. 
However, a good fit was also obtained by a single, neutral constant 
density reflector with twice the solar iron abundance. From those 
data it was unclear which interpretation was more plausible. 

Producing a physical model for the X-ray spectrum and spec- 
tral variability of MCG-6-30-15 is highly involved and there is no 
unique solution to this problem. In this paper we explore the impli- 
cations for models in which the spectrum is dominated by features 
from a photoionized ac cretion disc. Sp ecifically, we improve on the 
results of Ballantvne & Fabian 1 2001 ) and develop a model for the 
spectrum of MCG-6-30-15 in terms of reflection from two distinct 
sites on the disc and test this model against the observed spectral 
variability. An alternative model, in which the reflection spectrum 
is physically disconnect ed from the continuum d ue to relativistic 
effects, is discussed by Fabian & Vauahan (2003). We make use 
of data obtained from a 325 ks observation of MCG-6-30-15 by 
XMM-Newt on and fr om a simultaneous observation by BeppoSAX 
iFabian et all2002bl) . 

In the next section, the data and models are briefly introduced. 
In Section [3] we use the broad-band BeppoSAX data to constrain 
the reflection spectrum in a manner independent of the iron line. 
Section[4]presents the results of the ionized disc fits. The findings 
are discussed in Section|5|before conclusions are presented in Sec- 
tion[6] 



2 DATA AND MODELS 

Details of the observations are given in IFabian et alJ 12002b). For 
the present paper the data were processed using XMM-Newton Sci- 
ence Analysis System (SAS) v5.3.3, with source data extracted 
from circular regions of radius 35" from the EPIC MOS and pn. 
Background events were extracted from off-source regions. Events 
corresponding to patterns 0-12 (single-quadruple pixel events) 
were extracted from MOS and patterns 0-4 (singles and dou- 
bles) were used for the pn analysis, after checking for consis- 
tency with the data extracted using only single pixel events (pat- 



tern 0). Standard redistribution matrices (ml_r 6_all_15 . rmf ) 
for MOS1 (similarly for MOS2) and (epn_sw2 0.sdY9 . rmf) for 
the pn were used, and ancillary response files were generated using 
ARFGEN V 1 .48.8. 

The additional calibration information include d in this ver- 
sion of the SAS now allow the use of the pn data (in Fabi an et alJ 
(2002b) only the MOS data was analyzed). However, it was found 
there was a small but significant difference in photon-index (Ar ~ 
0.1) between the MOS and pn spectra even after removing the small 
amount of pile-up in the MOS data. Since the value of V obtained 
from the pn data was more consistent with that from the MECS 
and PDS instruments, we ig nored the MOS data in this analysis 
(see IFabian & Vaughanl 120031) for further discussion on the EPIC 
data). 

The con st ant d ensity ionized reflection models of 
iRoss & Fabianl 1 19931) and iRoss et alJ 1 19991) were used in 
the spectral fitting, as these are the most useful in parameterizing 
the shape of the spectrum. Models wh ere the accretion disc 
atmosphere is in hydrostatic balance (e.g., jNavakshin et alJl200ri 
IBallantvne et alJl200lt ; fRozariska et al]|2002l) were not used be- 
cause they are depend ent on the un known accretion disc structure. 
However, Ballantvne et al. 1200 ll) found that the hydrostatic 
models could be fit with dilute versions of the constant density 
ones. Therefore, there is not a significant difference in the shape of 
the reflection spectrum between the two cases. 

XSPEC vll.2 lArnaudl[T99d) was used for all spectral fitting, 
and we use the 90 per cent confidence level to compute errorbars 
for the fit parameters. 



3 THE REFLECTION STRENGTH 

As a first step toward determining the properties of the reflection 
spectrum in MCG-6-30-15, the high-energy BeppoSAX PDS data 
were examined in order to constrain the strength of the Compton re- 
flection "hump." In this analysis models of the reflected continuum 
were compared to the PDS data over the 14-200 keV range, which 
includes the peak of the reflection spectrum. The MECS data were 
included in the fitting to better constrain the slope of the under- 
lying continuum. Only MECS data in the ranges 2.5-3.5 keV and 
7.5-10.0 keV were used as these should be dominated by the pri- 
mary X-ray continuum; the spectrum is complicated by absorption 
below 2.5 keV and the strong, broad iron line in the ~ 3-7 keV 
range. The normalizat ion of the PDS data was t ied to 0.86 times 
the MECS value iFiore. Guainazzi & Grandill999l) . The EPIC data 
were not included in these fits as the much higher signal-to-noise 
ratio of these data would mean the model parameters would be de- 
termined almost entirely by the low-energy EPIC data, not the PDS. 

A simple power-law model for the broad-band continuum pro- 
vided a poor fit to the data (x 2 = 170.7/46 de grees of freedom, 
d.o.f.), with a photon index of V = 1.85. The residuals from this fit 
showed the PDS data were systematically in excess of this model in 
the 14-40 keV range, a strong indication of the presence of an ad- 
ditional Compton reflection continuum, which peaks at ~ 30 keV. 
Allowing for an exponential cut-off in the power-law continuum 
did not improve the fit (x 2 = 170.7/45 d.o.f.) and gave a best- 
fitting cut-off energy E cut > 500 keV, i.e. above the band-pass 
of the data, implying that the high energy spectral curvature is not 
the result of a roll-over in the power-law continuum. This fit is il- 
lustrated in Fig.Q The earlier BeppoSAX observation reported by 
Guain azzi et alJ |1999) also clearly detected the reflection contin- 
uum of MCG-6-30-15. 
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Figure 1. Top panel: BeppoSAX PDS and MECS data. Middle panel: resid- 
uals from a fit using an exponentially cut-off power-law model. Bottom 
panel: residuals allowing for reflection as described in Sectionl3l 

The broad-band spectrum was therefore modeled in terms 
a powe r-law plus Compton reflectio n continuum. The PEXRAV 
model ( Magdziarz & Zdziarski 1995) was used to compute the 
spectrum. Specifically, this computes the spectrum resulting from 
an exponentially cut-off power-law continuum incident on a slab 
of optically thick, neutral material, accounting for Compton reflec- 
tion and bound-free absorption. Unlike the lRoss & FabiaJ 1 19931) 
ionized disc model, PEXRAV does not account for discrete emis- 
sion features (specifically, fluorescence and recombination emis- 
sion). However, it does account for the angular dependence of the 
reflected spectrum and the reflection of incident photons with en- 
ergies ;> 100 keV (wh ich requires relativistic corrections to the 
scattering cross-section; White. Lightman & Zdziarski 1988). This 
model should therefore be appropriate for the energy ranges in- 
cluded in the fitting, which extend above 100 keV and ignore re- 
gions containing strong atomic emission features. 

Including a PEXRAV reflection component in the model im- 
proved the fit substantially (x 2 = 71.7/44 d.o.f.) and gave the 
following fit parameters: T = 2.17 ± 0.10, E cut > 143 keV and 
a normalization of the reflected spectrum relative to the primary 
continuum of R = 2.7i 'g, where R = Q./ 2n. The inclination 
angle of the reflector was fixed at i = 30° iTanakaet ailll995t 
Fabian et al. 2002b) and the abundances were assumed to be solar. 
This fit suggests rather strong reflection; an isotropic continuum 
source illuminating a flat accretion disc would be expected to yield 
a relative reflection strength of R ~ 1 (i.e. the reflector subtends a 
solid angle ~ 2n as seen by the continuum source). 

The exact strength and shape of the observed reflection spec- 
trum does however also depend on the abundances of the reflector, 
its inclination and and its ionization state. The effects of these pa- 
rameters on the measured value of R were investigated as follows. 
The PEXRAV model was re-fitted to the data with the Fe abundance 
fixed at 3x and then 7x solar and also with the Fe abundance as a 
free parameter. The 3 x solar model made the largest improvement 
to the fit (compared to the solar case) and the best-fitting abundance 
value was 3.5^}'g. This gave \ 2 — 63.6/44 d.o.f. and the reflec- 
tion fraction was constrained to be R > 2.6 (this fit is shown in 
the bottom panel of Fig.0. Allowing the abundances of elements 
lighter the Fe to vary did not improve the fit further. The Fe abun- 
dance was therefore fixed to be 3 x solar in the remaining PEXRAV 
models. iLee et all ll999l) previously presented evidence for an en- 
hanced iron abundance in MCG-6-30-15. 



The apparent strength of the reflected spectrum, for a given 
value of R, is also dependent upon the inclination angle of 
the reflector with respect t o the line-of-sight (e.g., Fig. 5 of 
Maedziarz & Zdziars ki 19951) . Allowing the inclination angle to be 
a free parameter gave only a small improvement in the quality of 
the fit (x 2 = 63.4/43 d.o.f.) with i < 80° and similar constraints 
on the reflection strength R > 2.8. 

The final possibility explored to explain the high value of R 
was ionization of the reflector. This was performed by using the 
PEXRIV model and allowing the ionization parameter £ to vary 
(with 3x solar Fe and i = 30°). This gave only a slight improve- 
ment in the fit (x 2 = 61.8/43 d.o.f.) with £ < 200 and R > 2.8 
(this £ has been corre cted for th e different energy ranges used in 
PEXRIV and the Ross & Fabian 1 1993) model). The low value of 
£ means that the lack of Comptonization of the Fe K edge in this 
model will not affect the limit on R. Thus, even allowing for the 
possible effects of non-solar abundances, angular-dependence and 
ionization on the reflection spectrum, the strength of the reflection 
measured using the high-energy PDS data appeared to be R > 2. 
This constraint on R is independent of the possible presence of a 
roll-over in the power-law continuum, since the cut-off energy of 
the power-law was a free parameter in the above fits. 



4 IONIZED REFLECTION RESULTS 

The BeppoSAX constraint on R is difficult to explain in terms 
of the standard disc-corona geometry which predicts R ~ 1. 
A large reflection strength can be generated if the X-ray emit- 
ting plasma is mo ving toward the accretion disc at mildly rel- 
ativistic velocities jRfivnold^^FabmrJh^97l: iBeloborodovl 1 1 9991 
Malzac, Beloborodov & Poutanen 2001). However, this model 
would predict an extremely soft V > 2.2 spectrum for MCG-6- 
30-15, which is not obse rved. Ano ther method of obtaining a large 
value of R was given by iFabian et alj {2002a ) who considered the 
resu lts of reflection from a lumpy or c orrugated accretion disc (see 
also iRoss. Fabian & Ballantvnd l2002). In this scenario, the opti- 
cally thick disc may subtend a solid angle larger than 2-n as seen 
from the X-ray source. Again, softer spectra are naturally predicted 
from such geometry, as the hard X-ray source 'sees' more UV pho- 
tons from the disc. 

When fitting the A .SCA spectra of MCG-6-30-15, 
iBallantvne & Fabian 1200 ll) found that it was very difficult 
for an ionized Fe Kq line to fully account for both the strength 
and width of the observed feature. On the other hand, an ionized 
reflection continuum provided a good fit to the obse rved one, in- 
cludin g the red wing of the Fe Ka line. Thus. lBallantvne &~F abian 
1200 ll) were forced to consider a two component model, where the 
red wing was fit by the ionized reflector, and the line core by a 
neutral diskline at 6.4 keV. Misra 1 2001) has also considered a two 
component model of the MCG-6-30-15 iron line. A superposition 
of more than one reflection event may also account for the observed 
reflection strength. In this section, we investigate the validity of 
such a model with our high signal-to-noise XMM-Newton data. 

4.1 The time-averaged spectrum 

MCG-6-30-15 has a strong warm absorber that significan tly af- 
fects the continuum below 2-3 keV j Njmdra & Pounds! Il992t 
iRevnolds etail 1 1995b ILee et all 1200 it ISako et alJ l2002t) . There- 
fore, the use of the full pn spectrum would necessitate a com- 
plete model of the warm absorber from analysis of the Reflection 
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Figure 2. The summed double reflector model fit to MCG-6-30-15 between 
2.5 and 80 keV. The top panel shows the total spectrum in EFg space (see 
Figurel5lfor a plot showing the model components), and the lower panel are 
the residuals in units of standard deviations. The pn data are shown with the 
open symbols, and the fit parameters are displayed in TableHI The residuals 
at 20-30 keV are a result of the reflection spectra being angle-averaged, 
which reduces the strength of the reflection hump I Magdziarz & Zdziarski 
1 995 ) for inclination angles < 50 degrees. 



Grating Spectrometer (RGS) data. However, the strongest reflec- 
tio n features are not affected b y t he warm absorber, so , follow- 
ing B allantyne & Fabianl l200ll) and lFabian et alj<2002bl) . we only 
considered data above 2.5 keV. The MECS and pn data were cut- 
off at 10 and 12 keV, respectively The P DS data were cut-off at 
80 keV because the Ross & Fabianl il9 93) models have a high en- 
ergy limit of 100 keV. All fits included the Galactic absorbing col- 
umn of 4.06 x 10 20 cm" 2 iElvis. Wilkes & Lockmarll989h . Data 
from the three instruments were fit simultaneously with the normal- 
izations free except for the PDS which was fixed to be 0.86 times 
the MECS value as in Sect.|3] 

The model consisted of the sum of two reflection spectra with 
solar iron abundance. The first of these components was relativis- 
tically blurred with a lLaorl ll99lh kernel that had a power-law 
emissivity (oc r - ' 3 ) with the index frozen at /3 = 3. The second 
component was also blurred, but with a flatter index of (3 — 2 ap- 
propriate for the case when the disc is warped and emission from 
large radii is non-negligible. Both reflectors were illuminated by 
the same power-law spectrum, but this primary emission was in- 
cluded only in the first component where it was added in so that R 
was unity. Thus, the second component was reflection dominated. 
In this way, a total observed R > 2 can be obtained, a s requir ed 
by the BeppoSAX data (Sect.|3). Following Fabi an et alJ l2002bl) a 
Gaussian with zero width was included at 6.9 keV to fit an addi- 
tional spectral feature. 

The best fit obtained with this model is shown in Figure|2|and 
TableQ We find a very good fit to the data (x 2 /d.o.f.=1600/1879) 
with an ionized reflector located at about 5 r g and a neutral re- 
flector outside of ~70 r„. T he photon index is consistent with the 
results of lFabian et all l2002bl) (T = 1.95; see their Table 1). The 
inner radius found her e (~ 5 r Q ) is farther out than that found by 
IWilms et alJ <200ll) or lFabian et ail |2002b) (~ 2 r g ). This is be- 
cause the ionized reflection model has both a curved spectrum and 
an intrinsically broadened Fe Ka line (due to Comptonization). 



The ionization parameter of the inner reflector (log£ = 3.82) is 
such that He-like iron is the dominant ion at the surface of the disc. 
Therefore, the rest energy of the Fe Ka line is at 6.7 keV, but is 
gravitationally redshifted down to below 6 keV. 

Examining the residuals to this fit in Fig.|2|shows that the re- 
flection hump is n ot accurately repr oduced. This problem is due 
to the fact that the R oss et all ^99) reflection models are angle- 
averaged rather than computed for a specific inclination angle (as is 
done in PEXRAV). For angles < 50 degrees, angle-averaging results 
in a w eaker reflection hump at 20-30 keV tMagdziarz & Zdziarskil 
Il995t). A possible secondary effect is the hard 100 keV cut-off in 
the R oss et alJ fl999) models causing the spectrum to turn down- 
wards faster than is observed. However, cutting off the PDS data at 
lower energies (say, 40-50 keV) makes little impact on the model 
parameters. Despite these issues, Fig. [2] illustrates that more than 
one reprocessing event can account for the strong reflection hump 
observed in the BeppoSAX data. 

Allowing [3 for the inner reflector to be a free parameter does 
not improve the fit for any value between 10 and — 10 because the 
inner and outer radii are so close together (see below). Doing the 
same for the outer component improves the fit only slightly with 
0.641 ^ /3 ut ^ 2.268. This result emphasizes that the model 
requires a flat emissivity at large radii, consistent with a strongly 
warped or concave disc. Allowing the reflection fraction to be free 
does drop the \ 2 by 8, which is significant at 99.8% according to 
the F-test. The best fit value of R in this case is 1.98^ ' 44. This 
additional reflection only affected the fit around the Fe Ka line, 
and did not alter the residuals at higher energies. 

In this model, the large red wing of the Fe line and the 2- 
10 keV continuum are due to the inner ionized reflector, while the 
blue core of the line and the reflection hump is dominated by the 
outer reprocessor. The most striking feature of this fit is that the 
radii over which the inner component originates is restricted to a 
narrow annulus only 0.08 r g wide. In Sect. l5.2l we discuss the im- 
plications of this model and attempt to evaluate if it is physically 
plausible. 

4.2 Time resolved spectra 

In order to investigate whether the double reflector model can 
explain the spectral variability properties of MCG-6-30-15, we 
split the XMM-Newton dataset into 32 10 ks segments (see Fig- 
ure Source and background spectra were extracted from the pn 
for each interval. Due to the ~71% observing efficiency of the 
pn IStriider et alj|200lh . the exposure time of each spectrum was 
-7.1 ks. 

Each of the 32 spectra were fit with the double reflection 
model between 2.5 and 12 keV using the time-averaged results (Ta- 
ble0 as a template. All parameters were allowed to float, except for 
the inner and outer radii of the distant reflector, the reflection frac- 
tion R, and the inclination angle i which were all fixed at the values 
listed in TableQ Keeping R fixed allows us to test whether the data 
are consistent with the innermost reflection component varying (in 
flux) along with the continuum (allowing for changes in the ioniza- 
tion parameter). The emissivity indices also remained unchanged 
from the time-averaged fitting. The model provided an acceptable 
fit (reduced \ 2 ~ 1) to the continuum and Fe Ka line for each 
spectrum. 

The results from this fitting are shown in Figure [4] The left 
panel shows the ionization parameter of the two reflectors as a func- 
tion of the total 2-10 keV flux calculated from the model. This plot 
clearly shows that the inner reprocessor at ~5 r g becomes more 
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Table 1. The parameters for the best fit time-averaged double reflection model shown in Figure l2l T is the photon-index of the incident power-law 
continuum, £ is the ionization parameter of the reflector, r; n is the inner radius of the reflector in r g , r ou t is the outer radius (also in r g ), and i is the 
inclination angle in degrees. The EWs of the 6.4 keV line and the narrow 6.9 keV line are given in eV. The EW of the broad Fe Ka line includes both 
components to the line profile. The inner reflector had R fixed at unity, and an emissivity of (3 = 3. The outer component was reflection dominated and 
had an emissivity of j3 = 2. p means the parameter pegged at its lower limit. 

Inner reflector Outer reflector 

T log£ r in r out logi; r in r out i EW 6 . 4kcV EW 6 . 9kcV x 2 / d -°-f- 



1.92±%°l 3.82+°;°^ 4.93±0.01 5.012+°;°°| p < 1.07 71+g 1822±^!f 31.6+% 450 29±f 1600/1879 

301 302 303 
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Figure 3. XMM-Newton EPIC pn light curve (0.2-10 keV) of MCG-6-30-15. The dotted lines denote the 10 ks segments discussed in the text. Source and 
background spectra were extracted from each segment. The total good exposure time of each spectrum was ~7.1 ks due to the live-time of the pn in small 
window mode. The count rates in this figure have not been corrected for this scale factor. 



ionized as MCG-6-30-15 becomes brighter. The relationship be- 
tween £ and the flux is roughly linear over most of the plot. The 
red wing of the Fe Ka line (which is fit by this reflector) will there- 
fore be weak at high flux levels. In contrast, the ionization param- 
eter of the outer reflector is uncorrelated with the 2-10 keV flux 
and remains small, consistent with predominantly neutral reflec- 
tion. The right panel in Fig.|4]plots T as a function of 2-10 keV 
flux. Here, we see that we have recovered the well known property 
that t he spectrum of MCG-6-30-15 softens as it becomes brighter 
(e.g JVaughan & Edelsorl200lLlShih et alfeOOA . 



with a broken power-law emissivity profile. In this model the ex- 
tended red wing of the line is produced by emission within 6 r B , 
while the core of the line profile (which was resolved) is produced 
from reflection beyond 6 r g and with a relatively flat emissivity 
profile. These analyses both suggest that the core of the line around 
6.4 keV is mostly produced by reflection off material at distances of 
~ 6 — 1000 r g , while the red wing of the line is produced by emis- 
sion from within 6 r g . Such an interpretation is consistent with the 
ionized reflection results found above. In this section we consider 
the two component ionized disc model in detail. 



5 DISCUSSION 

XMM-Newton has observed MCG-6-30-15 twice. The first obser- 
vation showed the broad iron line likely extends down below 4 keV, 
implying emissi on from within 6 r n , and also showed a narrower 
core to the line IWilms et all200lh . This apparent two component 
line may suggest that there is reflection from two distinct sites. 
The second reprocessing site is unlikely to be very distant (e.g., 
the molecular torus postulated in Seyfert unification schemes) as 
the high resolution Chandra HETGS spectrum showed little flux 
from an intrinsically narrow line component <Lee et all20o3) . The 
second, longer XMM-Newton observation also appea red to show 
line emission extending down to low energies. lFabianetalJl2002ri) 
modeled this in terms of neutral reflection from an accretion disc 



5.1 General comments 

In Section [4] it was found that ionized reflection models can fit 
the XMM-Newton data of MCG-6 -30-15 but only with some dif- 
ficulty. As was found previously iBallantvne & Fabianll200ll) . it 
is difficult for the predicted spectrum to account for the strength 
and width of the Fe Ka line, as well as the smooth continuum. 
The most tenable model required a secondary reflector from far- 
ther out on the disc, but, more worrying, the region of primary 
X-ray generation is confined to a very narrow annulus. Are there 
limitations in our reflection spectra which is resulting in these 
extreme models? Although it is difficult to answer this ques- 
tion definitively bec ause the de tailed physics of AGN is rela- 
tively unknown, the Ross & Fabianl ^93) reflection models have 
been very successful in fitting the X-ray spectra of many differ- 
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Figure 4. Left: Plot of the ionization parameter for the two reflectors versus the total 2-10 keV flux. The inner reprocessor clearly become more ionized as the 
source becomes brighter, while the distant reflector remains neutral. Right: Plot of the photon index (T) versus the total 2-10 keV flux. The source becomes softer 
as it brightens, which is in agreement with many other observations of MCG-6-30-15 using simpler models. 



ent Seyferts lBall antyne 1 Iwasawa & Fabian|200l|;[O rr et all200lk 
iDe Rosa et alJl2002t iDe Rosa. Fabian & Piral2002ll . On the com- 
putational side, the code has been compared with similar angle- 
averaged reflection calculations I Peauignot et al. 2001) and there 
was good agreement at hard X-rays, particularly around the Fe Ka 
line. The models used here assume reflection from a constant den- 
sit y slab, rather than a variabl e density atmosphere, as was done 
in Ballantvne & Fabian 1 2001). In these variable density models 
weaker Fe Ka lines are expected due to the smaller density at the 
surface of the disc; therefore, a constant density slab, although per- 
haps unrealistic, is the "best case scenario" for producing strong 
Fe Ka lines. Thus, it seems unlikely that this extreme fit is a result 
of limitations in the reflection models. 

It is worth po inti ng out that, consisten t with the results of 
IWilms et alJ 1200 ll) and lFabian elatl <2002b» . the Fe Ka line pro- 
file in this model suggests a concentration of emission from within 
~5 r g . It is unknown whether this is from the plunge region of a 
Schwarzschild black hole, or from the disc around a spinning Kerr 
hole. In either case, it seems as though additional physics must be 
included into our thinking of energy generation in AGN. 



5.2 Implications of a double reflection model 

5.2.1 Observational testable predictions 

In Figure [5] the double reflection model is broken into its compo- 
nents and plotted between 0.3 and 80 keV. Absorption from the 
Galaxy has been removed so that the soft X-ray spectral features 
can be seen unhindered. The total 2-10 keV flux given by the model 
is 4.36 x 10 -11 erg cm~ 2 s — l , of which only 6% is contributed 
by the distant component. When the energy range is extended to 
100 keV, the outer reflector provides 21% of the total flux. Thus, if 
this model is viable we would expect a decrease in rapid continuum 





ii 


















\ 1 1 n A- 1 ' \ 
jj'f ;j \ 






\ / 

s / 

I- / 

\ v / 


% Y 
\ 'i 

\ 

\ 










\- 




\ i 'J 




i, , I 






1 




10 





Energy (keV) 

Figure 5. The time-averaged reflection spectrum between 0.3 and 80 keV 
predicted by the double reprocessor model discussed in Sect. 14.11 Galactic 
absorption has been removed in order to show the soft X-ray lines predicted 
by the model. The solid line is the full model; the dashed line denotes the 
inner, ionized reflector; the incident T = 1.92 power-law is shown with 
the dotted line, and the outer, neutral reprocessor is the dot-dashed line. 
The 6.9 keV is also included in the plot. The noise at ~2 keV is a plotting 
artifact that did not affect the spectral fit. 

variability at higher energies, as the more distant emission would 
contribute to a larger fraction of the total spectrum. Furthermore, 
there may be a detectable lag in variations at high energies com- 
pared to those occurring between 2 and 10 keV. These predictions 
could be tested by Integral or ASTRO-E2. 

Another observationally testable aspect of this model are the 
soft X-ray spectral features. The ionization parameter of the outer 
reflector is constrained to be small enough for iron to be only 
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Table 2. Equivalent widths (EWs) of the soft X-ray emission lines predicted 
by the double reflection model when log £ = 1.0 for the outer reflector. 

Line EW (eV) 

SiKa 6 

Mg K« 4 

O VIII Ly« 9 

O VII 36 

N VI 6 

C VI Lyct 9 



weakly ionized, but in these cases recombination emission from 
ions of highly ion ized metals such as C, N and O are common 
at lower energies jRoss &Fabiar]|l993l) . As there is only an up- 
per limit on the £ of the outer component, we cannot predict the 
exact lines which may be present in the total spectrum. However, 
as an example, the EWs of the soft X-ray lines relative to the to- 
tal predicted continuum are shown in Table [2] for the case when 
l°g£ = 1-0- The strongest line is predicted to come from O VII 
with an EW of ~36 eV. Once an accurate model for the warm 
absorber has been determined it should be possible to search for 
such a line in the current dataset. The line is not extremely rel- 
ativistically blurred, so may even stand out in the RGS spectra 
and would be unrelated to the relativistic soft X-ray line model of 
iBranduardi-Ravmont et alJfcOOll) . 

In Sect. Pit was found that the ionization parameter of the 
inner reflector was correlated with the 2-10 keV flux and was al- 
ways greater than 1000. This implies that the strength of the iron 
line (which comprises the red wing of the obse r ved Fe Kg line) 
is ant icorrelated with the flux iMatt et ail 1 19961 : iBallantvne et alJ 
120021) . However, if the flux of MCG-6-30-15 drops to about one- 
half of the faintest level observed here so that log £ < 3, then the 
strength of the red wing correlates with £. This is another observa- 
tional testable feature. 



5.2.2 Geometry 

The geometry implied by the double reflector fit requires two dif- 
ferent changes in curvature in the accretion disc. First, the disc must 
become strongly concave at a radius of ~50-70 r g to provide the 
surface for the second reflection component. A standard Shakura- 
Sunyaev gas-pressure dominated accretion di sc is already mildly 
conca ve with the scale height H cx r 9//8 1 Shakura & Sunvaev 
Il973l) . but a warp much larger than this is re quired to substan - 
tially increase the reprocessed flux (see Fig. 2 in| Blackmail 1 999). 
A possible mechanism m ay be warping due to radiation pressure 
on the disc (Pringle 1996, 1997), however the typical transition ra- 
dius for this instability is ~ 1000 r a for an AGN disc JPringld 
I l997l). The Bardeen-Pe tterson effect feardeen & Pettersonlll97a 
Kumar & Pringlell985l) . which aligns the inner part of an accretion 
disc so that it is perpendicular to the spin axis of the central object, 
is another possible method to warp a disc if the outer regions are 
not in the equato rial plane. The transition radius for t his mecha- 
nism is < 100 r g llvanov & Illarionovll997LlNelson & PapaloizovJ 
2000), so the Bardeen-Petterson effect can account for the outer 
disc geometry implied by the double reflection model. Of course, a 
spinning black hole would be needed to cause this effect. 

Another constraint on the geometry is provided by the inner 
reflector. Here, the emission is confined to a narrow ring at ~5 r g 
(see Sect. 15.2.41 and there is no reprocessing between that point and 
the outer disc, ten times further away. This implies that the inner X- 




Figure 6. A schematic representation of a possible geometry implied by the 
double reflection model. X-rays are produced at the inner edge of the accre- 
tion flow which is thicker than the rest of the inner disc (possibly because of 
large radiation pressure support). Reflected emission from this environment 
is greatly blurred and highly ionized. The increased height shields the inner 
disc from much of the direct radiation, and it is not until the disc becomes 
more concave at large radii that it produces more reflection. 

ray source is hidden from most of the inner accretion disc, which 
requires the scale height to rise quickly at ~5 r g . A schematic dia- 
gram illustrating this geometry is shown in Figure|6| The increase 
in disc thickness could be due to an increase in the radiation pres- 
sure support in the innermost region of the accretion flow. 

5.2.3 Fe Ka variability 

The most intriguing aspect of the double reflector model is its po- 
tential to possibly explain the puzzling variability of the Fe Ka line 
profile llwasawa et all 19961 1999). Analysis of the long (4. 5 day) 
1994 ASCA observation of MCG-6-30-15 bv llwasawa etall<199a) 
considered the variability of the line core separately from the broad 
red wing. They found that on timescales greater than a few 10 4 s, 
the narrow component seemed to increase with the continuum flux, 
while the strength of the broad component decreased. On shorter 
timescales of a few 10 3 seconds the broad component varied, but 
the narrow core was consistent with staying constant. Reinforcing 
these trends was the spectral evidence provided by examining the 
data during flares or dips in the light curve. llwasawa etlril 1 1996) 
found that during a bright flare, the spectrum (integrated over 36 ks) 
showed a strong line core, but a very weak red wing. In contrast, a 
15 ks deep minimum at the end of the observation showed a strong 
red wing and a weak blue core. Another interesting line profile 
was found during a bright, ho ur long flare in t he 199 7 long ASCA 
observation of MCG-6-30-15. llwasawa et alj 1 19991) reported that 
the spectrum extracted from the flare had line emission only from 
within 5 r g . There was little evidence for a core at 6.4 keV, with a 
limit on the EW of 60 eV. 

A long 400 ks RXTE observation also illustrated the strange 
behaviour of the Fe Ka line in MCG-6-30-15. Although the PC A 
on RXTE cou ld not reso lve the line int o a core and wing, these data 
were used bv lLee et alJ 1 200$ and lRevnoldsl feOOOl) to search for 
variations in the line flux with the continuum. Both studies found 
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that much of the line flux was constant and unresponsive to large 
continuum chan ges down to timescales as low as perhaps 0.5 ks. 
Vau ghan & EdelsorJfcOOll fit these data on the orbital timescale of 
RXTE and found variations in the line flux, but these were also un- 
corrected with the continuum. Fabia n et alj l2002bl) found a simi- 
lar result using this XMM-Newton observation. These authors sub- 
tracted a 10 ks low flux state spectrum from a 10 ks high flux state 
one and found that the difference was well fit by a power-law. The 
constancy of the Fe Ka lin e flux, which has b een seen in other 
Seyfert Is (e.g., NGC 5548: IChiang et alJl2000h . is hard to under- 
stand if the entire line arises from close to the primary X-ray source. 

Analysis of time-resolved XMM-Newton spectra using the 
double reflection model found that the inner component became 
more ionized when MCG-6-30-15 was brighter (the left panel of 
FigHJ. Recall from Fig.|5|that the spectrum between 2.5 and 12 keV 
is dominated by emission from the inner zone. Also, the shape o f 
an ionized reflection spectrum changes with £ iRoss et all 1999). 
Thus, the change in ionization parameter found in the inner reflec- 
tor is due to variations in the shape of the spectrum, and not just the 
Fe Ka line. This is not necessarily an obvious result. MCG-6-30- 
15 softens as it becomes brighter (see the right panel of Fig.|4j, and 
a steeper illuminating power-law will decrease the ionization state 
of the reflector at a constant irradiating flux. Therefore, in order for 
£ to increase the incident flux must grow at a rate fast enough to 
offset the changes in the spectral slope. 

As the ionized reflector is nearly coincident with the X- 
ray source in this model, reverberation would be large on short 
timescales, adjusting the shape of the line. Therefore, assuming 
a constant R, the changes in the Fe Ka line profile observed by 
ASCA can be explaine d by changes in the ionization state of the in- 
ner reflector (see also Revnolds 2000). When the source is bright, 
the reprocessing zone becomes highly ionized and the amplitude 
of the red wing diminishes. On the other hand, when the source 
is in a low state, the ionization parameter drops an d there is more 
emissi on from the red wing explaining the results o fllwasawa et alJ 
(1996). However, can the blue core from the distant, neutral re- 
flector respond on the correct timescale to explain its variations? 
The best fit inner radius for the outer reflector is 71 r g (Table0. 
This corresponds to a light crossing time of 3.5M7 ks, where 
M7 = M/10 7 Mm. Ilwasawa et al] 1 19991) saw no blue core in 
their 1997 ASCA flare spectrum which lasted about an hour - the 
same timescale estimated above. Therefore, increasing the black 
hole mass by only a small a mount w ill comfortably explain this 
observation. Yet Fabi an et alJl200 2b) found that the difference be- 
tween a high-state and low-state spectrum of MCG-6-30-15 is con- 
sistent with a power-law. The conclusion from that test is that the 
Fe Kq line must not have changed strength or shape between the 
two states. 

To investigate the line variability predicted by the ionized disc 
model more closely, we show in Figure0three examples of a dif- 
ference spectrum obtained by subtracting the best-fitting double re- 
flector model for a low flux state segment from the model of a high- 
state segment. Dips in these spectra correspond to regions where 
the feature was stronger in the low flux state than in the high flux 
state, and bumps correspond to areas where the high flux model was 
stronger than the low flux model. The difference spectra obtained 
from the data were unfolded and plotted in Fig.Qto illustrate their 
statistical quality. 

The bottom panel shows the difference between the models 
for Rev. 301:e and 301:f (see Fig.|5J- At the end of 301:e there is 
sudden drop from a 2-10 keV flux of 5.21 x 10~ n erg cm -2 s _1 
to 3.07 x 10~ n erg cm -2 s _1 in ~ 1000 s. The ionization param- 
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Figure 7. Unfolded spectra obtained from subtracting the best-fitting double 
reflector model from a low flux state interval from the model of a high state 
segment, (a) The result when segment Rev. 303:j and 303:b (see Fig. l3l 
were used as the low and high flux states, respectively. These segments are 
separated by 70 ks. (b) The result when segment Rev. 302:c and 3()2:b were 
used as the low and high states. These intervals are adjacent to each other, 
(c) The result when segment Rev. 301:f and 301:e were used as the low and 
high states. These intervals are also adjacent to each other. See the text for 
the discussion of the interpretation. 
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eter of the inner reflector fell from log£ = 3.96 to 3.61 over that 
interval, as did the photon index: T = 1.98 to 1.87. The overall 
negative slope of the difference spectrum is due to that change in 
photon index. At the energy of the red wing of the Fe Ka line there 
is a slight deficit with respect to the continuum. This implies that 
the red wing was stronger in the low flux state spectrum (as a result 
of the lower £) than in the high flux state spectrum. This feature 
is weak and very broad and is difficult to detect in the data. The 
difference spectrum also shows a hump at the position of the line 
core. This means that the 6.4 keV core was weaker in the low flux 
state than in the high flux state. If the distant reflector had equili- 
brated at the high flux state in segment 301 :e, it would have not yet 
responded to the rapid decrease in flux unless the black hole mass 
was few x 10 6 M©; alternatively, it may be reacting to changes 
before the flare in 301 :e, as it lasted only 10 ks (see Fig.[3}. 

The middle panel of Fig.0plots the difference between Rev. 
302:d and 302:c, which illustrates the case of the source moving 
from a "steady" low-state (Rev. 302:a,b,c all have roughly the same 
flux) to a bright flare that occurs at the end of the Rev. 302:d. Fig.|5| 
shows that the pn count rate increased by nearly a factor of 2 over 
the 10 ks interval. The difference spectrum of these two segments 
shows a large (factor of 10) drop at ~ 6.5 keV implying that the 
core of the Fe Ka line was stronger in 302:c than in 302:d, opposite 
to what is expected if the line core was responding to the continuum 
on timescales of ~ 10 ks. 

The top panel of Fig. shows the difference spectrum be- 
tween two segments which are separated by 70 ks. The ioniza- 
tion parameter for the data of Rev. 303:b is log£ = 3.95 while 
it is log£ = 3.57 for segment 303:j. Thus, the red wing of the 
Fe Ka line, which originates from the ionized reflector should be 
stronger in the low flux state spectrum. This does seem to be the 
case, as the difference spectrum shows the profile of the red wing 
as a deficit between 3 and 6.4 keV which means that this part of the 
spectrum is stronger in the low flux state spectrum than in the high 
flux state. The core of the line is marginally weaker in the low state 
than in the high flux state, consistent with the lower flux level. Yet, 
both of these features are <10% from the continuum, and again, 
as seen from the data points, difficult to detect unambiguously. The 
2-sigma upper-limit on the EW of any Gaussian feature at 6.5 keV 
in the data difference spectrum of 303 :b and 303 :j is 120 eV (it is 
230 eV for 301 :e and 301 :f). 

The conclusion from these tests is that this model does pre- 
dict changes in the line profile on short timescales, in particular 
around the red wing. However, these changes will be slight due to 
the breadth of the feature and that is why they have not been de- 
tected before. Unfortunately, higher signal-to-noise data at small 
timescales will be needed to check these predictions. Finally, the 
model also finds that small changes in the 6.4 keV line core are 
consistent with the data, but we are unable to determine a timescale 
for the chan ges. Such variations have not been seen before (e.g., 
iRevnold J20 00 ). and will need to be found for this multiple reflec- 
tion model to be viable, although the signature may be complex. 
For example, a warped accretion disc is subject to Lense-Thirring 
precession, and this would cause variability in the outer reflector. 
Although a numerical simulatio n is needed to fully m odel the dy- 
namics of such systems lArmitag e & NatarauuT|[l999| ). the p reces- 
sion period can be estimated from jBard^eri&Prttersorll975l) 




where a is the black hole spin parameter. Therefore, precession can 
only cause long timescale variability in this model. 



5.2.4 The inner annulus 

Thus far we have ignored the problem that this model predicts that 
the primary X-ray emission originates from within a tiny annulus 
at 5 r g . The fit converged to this value by fitting a notch in the 
spectrum at ~3.6 keV; in order for this to remain unblurred, the 
outer radius was forced to be close to the inner radius, which is 
constrained by the width of the red wing. Replacing the pn data 
with the MOS data does not affect this result. It is possible that 
this feature is an emission line from hydrogenic Argon from the 
outer reflector which is not included in the reflection models. Fur- 
ther work needs to be done to eliminate this possibility. 

The requirement for the emission to be constrained to an annu- 
lus is very extreme as viscous dissipation throughout the accretion 
disc should go as r~ 3 1 Sh akura & Sunvaevlll973l) . If most of the 
accretion energy is transported to a magnetic corona above the disc 
then it could only be efficient in this small annulus which seems 
highly unlikely given that magnetic fields will be present through- 
out the flow. In the 1997 flare spectrum reported by I wasawa et alJ 
1 1999), they also found that the line emission seemed to arise from 
an annulus around 5 r g . One explanation considered by these au- 
thors was that this emission was caused by a flare rotating along the 
very inner edge of the disc. But, as these authors pointed out, this 
implied a black hole mass greater than 10 8 Mq as well as theoreti- 
cal problems about the height of the flare. In the model considered 
here, emission from an annulus at 5 r g is required at nearly every 
point in the observation, and cannot be due to a transient event such 
as a flare. 

One possible explanation for this energy generation is a con- 



ncction between a spinning black hole and the accretion disc (e.g. 


Blandford & Znaiek 1977; 


Li 2000, 2002a b). This was considered 


by Wilms et alJ 


1 2001) to explain the steep emissivity profile found 



in an earlier XMM-Newton observation of MCG-6-30-15. Since 
emission is required from within 6 r g it is likely that the black hole 
is spinning, and is therefore a source of energy that can be trans- 
formed into radiation. Returning radiation is another method to in- 
crease the emissivity of the inn er disc iMartocchia. Matt & Karas) 
l2002tlFabian & Vaugharl20ol . 



6 CONCLUSIONS 

Making progress in the study of accretion discs and AGN can only 
be made by comparing models with data. Here, we have us ed the 
sophisticated ionized disc models of iRoss & Fabianl 1 19931) to fit 
the latest XMM-Newton and BeppoSAX spectrum of MCG-6-30- 
15. An ionized disc has often been considered as a possible expla- 
nation for the strange variability behavior of the broad Fe Ka line. 

We found that an ionized disc model can fit the data and ac- 
count for the spectral variability if there are two distinct reflection 
regions on the disc. This model requires a very specific geometry 
in order to be viable, but can be tested by searching for soft X-ray 
emission lines that are broadened to a few 10s of eV, but are not 
asymmetric, and, in the future, by the continuum variability prop- 
erties above 30 keV. An apparent problem of this scenario is that 
it constrains ~80% of the total X-ray flux of MCG-6-30-15 to be 
produced within a very narrow annulus only 5 r g from the black 
hole. This is a stringent requirement that is difficult to explain with 
standard AGN theory. Nevertheless, the advantage of the double 
reflection model is that it can naturally explain the strange Fe Ka 
line variability that has been previously observed. Time-resolved 
spectral fitting found that the inner reflector becomes more ionized 
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as the source brightens, weakening the red wing of the iron line. 
Changes in the 6.4 keV line core are also found in the time-resolved 
fitting, but they do not necessarily correlate with the continuum. 
The model can also account for the R > 2 constraint implied by 
the BeppoSAX data by combining the reflection humps of multiple 
reprocessors. 

One similarity that t he model presented h ere has with previous 
work (Wilms et al. 2001; Fabian et al. 2002b) is the need for extra 
production of radiation within 6 r g . Therefore, irrespective of the 
number of reflectors or the ionization state of the disk, it seems 
likely there is some new physical process near the black hole of 
MCG-6-30-15 that needs to be untangled. 
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